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Abstract. We present a proper motion survey over about 200 square 
degrees towards the NGP, based on the material used for the construction 
of the GSC-II, that we are using to study the vertical structure and 
kinematics of the Galaxy. In particular, we measured the rotation velocity 
of the halo up to 10 kpc above the galactic plane traced by a sample of 
RR Lyrae and BHB giants for which radial velocities were used to recover 
the complete distribution of the spatial velocities. Finally, the impact of 
astrometric and spectroscopic GAIA observation are discussed. 



1. Introduction 



It is generally accepted that the Galaxy is constituted by four discrete main 
components, the bulge, the thin disk, the thick disk and the /la/o, which are 
characterized by distinctive stellar populations in terms of spatial distribution, 
kinematics properties, metallicity and age. A detailed knowledge of such galactic 
components is essential to achieve a complete description of the Milky Way, as 
well as of the various processes and evolutionary phases which occurred during 
the history of our Galaxy (and other galaxies too) and that are responsible for 
the existence and the properties of those components we observe today. 
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Ground based surveys providing photometry, proper motions and/or spec- 
troscopic observations have been carried out in the past to study the structure 
and kinematics of Galactic populations, and these will continue and be extended 
in the next years thanks to availability of all-sky catalogs such as GSC-II and 
USNO-B, based on multi-epoch photographic surveys, as well as other current 
and future photometric and spectroscopic surveys (eg. 2MASS, DENIS, SDSS, 
EIS, HES, HK, VST, etc.) that benefit from the availability of dedicated scan- 
ning and large field cameras as well as large multi-fiber spectrographs (eg. 2dF, 
6dF, FLAMES). 

Clearly, the GAIA mission will provide an enormous contribution to the 
understanding of the formation and evolution of the Galaxy, thanks to very ac- 
curate astrometric parameters complemented by photometric and spectroscopic 
measurements which will permit the direct determination of the 6D phase space 
distribution and the chemical abundance of large and complete samples of stellar 
tracers belonging to the different galactic components. 

Here we describe a new project which combines astrometric, photometric 
and spectroscopic data in order to investigate the kinematics of the outer halo 
by means of velocities derived from proper motions and radial velocities for a 
sample of RR-Lyrae and BHB giants towards the NGP. 

2. NGP survey: proper motions and radial velocities 

At the moment we have surveyed about 200 square degrees towards the NGP, 
and produced positions, proper motions and photographic photometry for about 
500 000 objects down to plate limits (Rp < 20.5). 



Table 1. Plate material 



Survey 


Epoch 


Pixel 


Band 


Emulsion + Filter 


POSS-I E 


1950-1956 


25 


yum 


E 


103a-E -I- red plexiglass 


POSS-I 


1950-1955 


25 


/im 





103a-O unfiltered 


Quick V 


1982-1983 


25 


ium 


Vi2 


IlaD-l-Wratten 12 


POSS-II J 


1988-1996 


15 


/im 


Bj 


IllaJ + GG385 


POSS-II F 


1989-1996 


15 


//m 


Rf 


IIIaF + RG610 


POSS-II N 


1990-1998 


15 


fim 


In 


IV-N + RG9 



Our material consists of 6.4° x 6.4° Schmidt plates from the Northern pho- 
tographic surveys (POSS-I, Quick V and POSS-II) carried out at the Palomar 
Observatory (Table 1). All plates were digitized at STScI utilizing modified 
PDS-type scanning machines with 25 fim. square pixels (1.7 "/pixel) for the first 
epoch plates, and 15 lira pixels (1 "/pixel) for the second epoch plates. The digi- 
tal copies of the plates were initially analyzed by means of the standard software 
pipeline used for the construction of the GSC-II (see e.g. Lasker et al. 1995 or 
McLean ct al. 2000). The pipeline performs object detection and computes pa- 
rameters and features for each identified object. Further, the software provides 
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Figure 1. Left panel. Formal errors of fitted proper motions as a 
function of the magnitude (all stellar objects). Right panel. Vector 
point diagram (VPD) of a sample of 15 QSO's, with 1 a error bars. 
Weighted means are (/XaCOs5) = —0.02 ± 0.23 mas/yr and (us) = 
+0.33 zb 0.28 mas/yr, respectively. Both plots refer to the POSS-II 
field no. 442. 



classification, position, and magnitude for each object by means of astrometric 
and photometric calibrations which utilized the Tycho2 (H0g et al. 2000) and 
the GSPC-2 (Bucciarelli et al. 2001) as reference catalogs. Accuracies better 
than 0.1-0.2 arcsec in position and 0.15-0.2 mag in photographic magnitude are 
generally attained. Relative proper motions were derived by applying the pro- 
cedure described in Spagna et al. (1996) and afterwards transformed to the 
absolute reference frame forcing the extended extragalactic sources to have null 
tangential motion. As shown in Figure 1, the typical precision (cr^ ~ 3 mas/yr 
down to Rf — 18) has been estimated from the formal errors of the fitted proper 
motions, while the zero point accuracy of the absolute proper motions have been 
tested by checking the mean motion of a set of known QSO's that give values 
smaller than 1 mas/yr on each component. 

Radial velocities and chemical abundances of the sample of RR Lyrai and 
BHB giants were derived by means of spectroscopic observations carried out with 
the 4m Mayall telescope at Kitt Peak and with the 3.5m TNG on La Palma. 
The data were processed with standard procedures and routines (IRAF), and 
typical errors are ctrv < 40 km s~^ and cr^pe/u] ~ 0-2 dex. 



4 



Spagna et al. 



NGP - CMD 



NGP - 17<Rp<18 




1 2 
Bj-Rp (mag) 



0.05 



0.00 



i -0.05 



-0.10 

-0.10 -0.05 0.00 0.05 
M„ cos 6 i"/yr) 



Figure 2. Color magnitude diagram {left panel) and vector point di- 
agram (right panel) for stars with 17 < Rp < 18. Both plots are based 
on data from the POSS-II field no. 442. 



3. The vertical structure 

The color magnitude diagram and the vector point diagram observed in one field 
of our NGP survey are shown in Figure 2. The observed distributions are the 
result of the complex mixture of the stars belonging to the various populations 
which are present towards high galactic latitudes. In particular, these may 
include: 

1. the flat and rapidly rotating old thin disk, 

2. the extended thick disk, including its metal weak tail, 

3. a flattened and slowly rotating inner halo, 

4. a spheroidal non-rotating outer halo, 

5. satellite debris and kinematics substructures. 

Actually, the physical properties of these components are not completely 
established and various problems arc still controversial. For instance: (a) the 
density scale factor, rotation and metallicity distribution of the thick disk; (b) the 
nature of the metal weak thick disk (MWTD) and its relation with the standard 
thick disk (satellite debris or initial phase of a dissipative formation?); (c) the 
halo velocity ellipsoid, the spatial distribution and axial ratio of the halo as a 
function of the distance, (d) the search of halo streams; (/) the determination 
of the luminosity and mass function of the faintest Pop.II stars (white dwarfs, 
late M dwarfs and subdwarfs). 

Due to the lack of trigonometric parallaxes for our and similar surveys, 
accurate distances and tangential velocities cannot be directly determined. In 
such a case it is convenient to analyze these large data-sets by comparing the 
observations against ad hoc Galaxy models which describe both the density and 
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Figure 3. Starcounts derived from plate XP444 (histogram) and com- 
pared with the distribution predicted by the Galaxy model {thick solid 
line), which includes thin disk (solid line), thick disk (dashed line) and 
halo (dotted line) components. 

kinematics, as for instance the Besangon model (Robin &; Oblak 1987), the IAS 
Galaxy model (Bachall, Casertano Sz Ratnatunga 1987), or the models developed 
by Mendez & van Altena (1996) and Chen (1997). To this regard, in Figure 3 
starcounts are compared against the Mendez's Galaxy model which has been 
extended to the photographic Bj and R^ magnitudes by Spagna (2001). 

Alternatively, using tracers with known brightness, it is possible to derive 
distance and space velocity by means of photometry, proper motions and spec- 
troscopic radial velocities. This approach has been adopted to investigate the 
kinematics of the outer halo traced by a set of 31 RR Lyrae and 65 BHB giants, 
as will be discussed in the following section. These objects are distributed be- 
tween 2 ^ Z ^ 12 kpc (V=12-16.5 mag) and accurate photometric parallaxes 
((Td/d ^ 10%) have been computed by means of the My vs. B-V relation for 
BHB giants derived by Preston et al. (1991), while for RR Lyrse stars we adopted 
the My as a function of metallicity or of Fourier components, as follows: 

My = 0.23 [Fe/H] + 0.92 (1) 
My = -1.876 log P - 1.158^1 + 0.821^3 + 0.448 (2) 

where Eq. 1 is from Cacciari (2002) and Chaboyer (1999), while Eq. 2 is based 
on the relation from Kovacs &; Walker (2001) where P is the period (in days). 
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Ai and A3 are the Fourier amplitudes (mag) of the fundamental and second 
harmonic components, respectively. The zero point has been calibrated by Kin- 
man (2002) with respect to the absolute magnitude (My = 0.6llo;J5) of R-R- 
Lyr derived from the HST/FGS parallax measured by Benedict et al. (2002). 
Finally, the extinction E(B-V) has been estimated from the maps of Schlegel et 
al. (1998). 



4. Halo rotation 

A retrograde rotation of the outer halo has been suggested by Majewski (1992) 
who measured a mean velocity {V) = —275 km s~^, which corresponds to a 
galactocentric retrograde velocity v^ot — —55 km s^^ adopting Vlsr = 220 
km s~^, from the analysis of a pure sample of halo subdwarfs at Z > 5.5 kpc 
towards the NGP. As shown in Table 2, this parameter is still controversial, in 
fact Carney (1999), after correcting the kinematics bias of his kinematically- 
selected subdwarf sample, found a net prograde rotation of about {V) = —196 
km s~^. On the contrary, Chiba & Beers (2000) measured a prograde rotating 
inner halo, Wj-ot — 20-60 km s~^, up to about 1 kpc, with a decreasing vertical 
gradient of dV/d\Z\ = —52 it 6 km s~^ kpc~^, while they did not detect any 
significant rotation above Z 1.2 kpc for very low abundance stars (—2.4 < 
[Fe/H] < —1.9), where contamination of thick disk stars should be negligible. In 
addition, their halo sample at larger distances (212 stars with 4 < Zyiax < 20 kpc 
and [Fe/H]< —1.5) still does not support any significant rotation: v^ot — ± 8 
km s~^. 



Table 2. Recent measurements of the halo rotation. 



Tracers 


N.ro Objects 


{V) (km/s) 


Reference 






Inner Halo 




RR Lyrae 


162 {\Z\ < 2 kpc) 


-210 ± 12* 


Layden et al. (1996) 


RR Lyra3 


84 {\Z\ < 2 kpc) 


-219 ± 10* 


Martin & Morrison (1998) 


RR Lyrae 


147 {\Z\ < 2 kpc) 


-217 ± 13* 


Gould & Popowski (1998) 


Subdwarfs 


{\Z\ < 1 kpc) 


-(160^200)** 


Chiba & Beers (2000) 




Outer Halo 




Subdwarfs 


21 (Z > 4 kpc) 


-275 ± 16 


Majewski (1992) 


Subdwarfs 


30 (^Max > 5 kpc) 


-196 ± 13 


Carney (1999) 


Subdwarfs 


212 (ZMax > 4 kpc) 


-220 ±8 


Chiba & Beers (2000) 


RR & BHB 


53 (2 < Z < 12 kpc 


) -285 ±17*'*** 


This survey 



(*) Heliocentric velocities. (Wrt. the LSR in the other cases.) 

(**) As a function of distance with a gradient dV/d\Z\ = —52 it 6 km s"^ kpc' 

(-2.4 <[Fe/H]< -1.9). 
(***) Preliminary results based on 18 RR Lyrae stars and 35 BHB giants). 
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However, reeently Gilmore et al. (2002), who carried out a spectroscopic 
survey at intermediate galactic latitudes of about 2000 F/G stars, revealed a 
significant excess of retrograde halo stars in their faintest magnitude bin (18 < 
V < 19.5) corresponding to a vertical distance \Z\ ~ 5 kpc. 

The fact that the velocities of halo stars do not match an exact gaussian 
distribution is well known (sec eg. Martin & Morrison 1998). How much this 
depends on the properties of the whole population or on the effects of kinematic 
substructures, such as satellite debris of ancient accretion events (eg. Helmi et 
al. 1999), remains to be established. 

As shown in Table 2, a preliminary analysis of our sample of RR Lyrae and 
BHB giants (Kinman et al. 2002) seems to support a retrograde rotation of the 
outer halo. In fact, we measured a heliocentric velocity {V) = — 285±17 km s~^, 
which corresponds to Vrot — —60 km adopting a solar motion with respect to 
the LSR, Vq = +5.25±0.62 km s'^, from Dchncn & Binncy (1998) and assummg 
Vlsr = 220 km s~^. This result is confirmed also by the separate analysis of the 
RR Lyrae and BHB giants, which both provide a retrograde rotation. 

How much this value may be affected by a velocity bias can be estimated 
by the level of systematic errors on proper motions which in practice give the 
main contribution to the U and V galactic components along this line of sight 
towards b ~ 90°. As reported in Sect. 2, we found systematic errors < 1 
mas/yr, from which velocity bias up to 20-30 km s~^ can be expected for such 
stars located at Z 5-6 kpc, on average. Clearly, this is a critical point which 
need further analysis. 

5. Kinematics simulation 

A Montecarlo simulation has been developed in order to compare the GAIA 
capability to recover the halo kinematics with respect to the velocity precision 
that can be attained by the current and future ground based surveys. For sim- 
plicity we considered an uniform spatial distribution with Pop.II-like kinematics, 
{crui o'y ) f^w) = (150, 100, 100) km s^^ and an observer rotating with a velocity 
of 220 km s~^ at a distance of 8 kpc from the galactic center. Implicitly we 
assumed to measure bright tracers, such as RR Lyrae, BHB and red giants, hav- 
ing apparent magnitude ~ 16 mag at a distance of about d ~ 10 kpc. The 
following cases were tested: 

• Case A (GAIA), with a-j^ = 10 ^as, = 10 ^as yr^^ (per component) 
and o"vr = 10 km s~^, from astrometric and spectroscopic observations; 

• Case B (ground-based surveys), with am-M = 0.2 mag, = 1 mas yr~-^ 
(per component) and civr = 10 km s~^, from astrometric, photometric and 
spectroscopic observations. 

Note that the two cases differ essentially in a factor 100 on the proper motion 
accuracy. In fact, at ^ 10 kpc the distance accuracy is the same for both cases. 
This points out the fact that, although radial velocities and distance moduli 
derived from ground based spectro-photometric surveys can attain similar ac- 
curacy, GAIA astrometry will Tiniquely provide accurate and reliable tangential 
velocities up to large distances. 
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Figure 4 shows the results of Montecarlo simulations with 5000 stars towards 
b = 90°. Top panels present the velocity errors (i.e. the differences between the 
observed velocity and the true value) for Case A, and those at the bottom for 
Case B. As expected, errors increase linearly with distance, and for GAIA the 
rms of AU and AV residuals vary from ^ 10 km s~^ at 5 kpc to 20-30 km s~^ 
at 10 kpc. For Case B the errors on U and V are about a factor 2-3 larger, while 
CTAW = o'Vr = 10 km s~^ in both the cases. Note that in Case A the error on 
the tangential velocity is in practice dominated by the distance uncertainty, a.^, 
while both the errors on proper motions and photometric parallaxes contribute 
to the velocity errors of Case B. Similar results are derived in other directions. 

These velocity errors should be compared with the typical motion of the 
halo stars (100-200 km s~^) for which GAIA will provide individual 3D spatial 
velocities with a significant signal-to-noise, (Ty/v, up to d 10-15 kpc, a dis- 
tance where ground based surveys are not able to measure reliable tangential 
velocities. GAIA will measure direct distances and velocities for large samples 
of halo tracers, selected in situ without kinematics nor metallicity bias, from 
which it will possible to determine accurately the halo velocity ellipsoid and its 
orientation. In particular, the GAIA //-arcsec level accuracy will be essential in 
order to take advantage of the y/N statistical factor and avoid the risk of ve- 
locity biases due to the presence of systematics errors affecting proper motions 
and parallaxes as discussed in Sect. 4. 

Finally, accuracy of the order of 15-20 km s"-*^, such as that attained by 
GAIA for d ^ 10 kpc, is the value requested to resolve kinematics substructures 
in the halo, as the satellite debris predicted by the hierarchical scenarios of 
galaxy formation. To this regards, Helmi (2002) estimated that GAIA will be 
able to recover 2/3 of the accretion events with a velocity accuracy better than 
20 km s~^ and aa/d < 20% down to 1^ ~ 18.5, or 1/2 of the events down to 
y ~ 15 mag. 

6. Summary 

Large field proper motion surveys (o"^ ~ 1-10 mas/yr) based on photographic 
surveys digitized with fast measuring machines (APS, MAMA, PDS, PMM, 
SuperCOSMOS) arc still useful tools for the study of the structure and kine- 
matics of the galactic stellar populations, especially when combined with radial 
velocities and chemical abundances from spectro-photometric observations. In 
particular, we have shown how bright halo tracers, such as BHB giants and RR 
Lyrae stars with V ^ 16-17, can be used to investigate in situ the halo kinematics 
up to z « 5-6 kpc via photometric parallaxes if proper motions with ~ 1- 
2 mas/yr accuracy and radial velocities with precision of about avr ~ 10-20 
km s~^ are available. For statistical analysis of large samples, systematic errors 
on the proper motions (and My too) are critical and can result in biased mean 
velocities and dispersions. These problems will be dramatically reduced by the 
astrometric parameters and spectroscopic Vr provided by GAIA, which will per- 
mit to determine (a) direct distances from parallaxes, (5) radial and tangential 
velocities with an accuracy of a few tens of km s~^ to d ~ 10-15 kpc for (c) un- 
biased and larger sets of tracers identified by means of its spectro-photometric 
and astrometric data. 
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Figure 4. Kinematics simulation towards NGP. Velocity residuals for 
Case A {top panels) and Case B {bottom panels). Solid lines show the 
precision level (ilcr) as function of the distance up to Z=15 kpc. Above 
each plot, the rms of the residuals of the whole sample is reported. 



